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Abstract 

We investigate some generalizations of the Hofstadter problem to higher dimensions 
with Abelian and non-Abelian gauge held configurations. We numerically show the hier- 
archical structure in the energy spectra with several lattice models. It is also pointed out 
the equivalence between the 7r-fiux state and the staggered formalism of Dirac fermion. 
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1 Introduction 

The fractal structure of the energy spectrum for the two-dimensional magnetic lattice system, 
known as Hofstadter's butterfly pQ, is one of the most exotic consequences of the quantum 
property of the low-dimensionality Such a fractal nature of the magnetic system can be ob- 
served in a simple tight-binding model in the presence of the magnetic field. This Hofstadter 
problem and its variants have been extensively discussed in various areas of physics and also 
mathematics. 

In this paper we extend the Hofstadter problem, which is originally considered in two 
dimensions, to higher dimensions with not only Abelian, but also non-Abelian gauge field 
configuration. So far some attempts to generalize it to the three-dimensional magnetic 
system were proposed in [H [31 HI E]. But its generalization to much higher dimensional 
system has been not yet studied in the literature. Actually, when we analyse topological 
matters in the lower dimensional system, the four dimensional point of view can be quite 
useful: topological insulators/superconductors in two and three dimensions 0E] are deeply 
connected to the four dimensional QHE [8 J through the dimensional reduction procedure [9,, 
I10| . The formalism discussed in this paper can be applied to arbitrary even dimensional 
lattice system in the presence of the magnetic field. Based on this formalism, we prove the 
7r-flux state |lll I12j . which has a gapless excitation in general, is essentially equivalent to 
Dirac fermion in arbitrary dimensions. Furthermore, when we apply the non-Abelian gauge 
field, we can obtain a hierarchical structure in the energy spectrum of the lattice model even 
in higher dimensions. 

This paper is organized as follows. In Sec. [2] we introduce arbitrary even dimensional lat- 
tice models with Abelian gauge field background configuration. We discuss the corresponding 
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Schrodinger equation to the magnetic system, and obtain the Harper's equation. We then 
comment on its connection to the non-commutative torus. We also provide a general proof 
for the gapless spectrum of the 7r-nux state, by referring to its equivalence to the naively 
discretized lattice Dirac fermion. In Sec. [3] we then consider the Hofstadter problem for 
the non-Abelian gauge field configuration. In particular the four-dimensional SU(2) theory 
is investigated as a fundamental example of the non-Abelian gauge theory. We show some 
numerical results of the model, and discuss the effect of the inhomogeneity of the background 
flux. Section [4] is devoted to a summary and discussion. 



2 Abelian gauge field models 

First generalization of Hofstadter problem is formulated in arbitrary even dimensions in 
the presence of U(l) gauge field. Before introducing a lattice model, we now consider the 
following background U(l) field configuration for the continuum theory [13] . 

^2s-i,2s = w s , F^y = for otherwise, (s = 1, • • • , r) . (2.1) 

We then apply the higher dimensional version of Landau gauge to this configuration, 

A 2s -i(x) = -lo s x 2s , A 2s (x) = 0. (2.2) 

Here the field strengths w s are quantized, 

2vr . , 

u s = j—n s , n s G Z, (2.3) 

L>2s~\l-<2s 

and thus the topological number is given by 



Q 



[ d 2r xe^-^ TxF ^ 2 ... F ^_ i4i2r =f[n s . (2.4) 

' J s=l 



This is regarded as the r-th Chern number. 

To discuss the Hofstadter spectrum, we then realize these configurations on the lattice. 



The gauge potential (2.2) is implemented by introducing the link variable, which can be 
regarded as the Wilson line, 

172,-1 (s) = e-*"' 3 *', U 2s (x) = 1. (2.5) 



Precisely speaking, we have to assign appropriate boundary conditions even for U 2s {x) |13j . 
In this case the field strength is restricted to the interval (0, 2ir) due to the lattice discretiza- 
tion |14] 115]. Furthermore they are characterized by the following fractions, 

^ = - H \ , (2.6) 
Is L 2s -iL 2s 

Here Lj is related to the system size as Li = La where a is the lattice spacing. Note that 
they satisfy < Pi/q% < 1. This fraction plays an essential role in the interesting spectrum 
of the model we discuss below. 
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2.1 Tight-binding model 

The lattice Hamiltonian with this background configuration is defined as 



(2.7) 



This is just the tight-binding hopping model, which describes the non-relativistic particle 
with the background magnetic field. Introducing the state = ^2 x ipx c t\0), we obtain the 



corresponding Schrodinger equation from the Hamiltonian (2.7) written in a second quantized 
form, 



2r 

5^ ^ t (x)V>x+/i + Ufa - fi)ip x . 



n=i 



(2.8 



We then solve the equation (2.8) by taking Fourier transformation. Remark the trans- 



lation symmetry of this model is slightly modified from the usual lattice model due to the 
background field, 

X 2s -1 ~ X2s-1 + 1, X2.s~X2s + q S - (2.9) 

Thus, writing the coordinate as X2 S = Q s y s + z s , the wavefunction is Fourier transformed as 



1 

V 



exp 



8=1 



2s -ix 2s -i + qsfasUs 



(2.10) 



where V stands for the effective volume of the system, V = L\ ■ ■ ■ L2 r /(qi ■ ■ ■ q r ). The 



Schrodinger equation (2.8) is rewritten in this basis as 

r 

$z+S + ^z-s + 2 COS ( fc 2s-l - U S Z S ) ljj z 



s=l 



(2.11) 



This is the higher dimensional version of the Harper's equation |16| . While an one-dimensional 
equation is obtained from the magnetic two dimensional system, we have a r-dimensional 
equation from the d = 2r theory. Furthermore the matrix size of this higher dimensional 
Harper's equation is N q x N q where N q = qi ■ ■ ■ q r . This means that the number of energy 
bands is just given by N q , and there possibly exist N q — 1 energy gaps. To discuss these 
energy gaps, one has to investigate a transfer matrix and its spectral curve associated with 



the Harper's equation (2.11). In this case, however, it is impossible to show these energy 



bands are totally gapped, due to its high dimensionality: there are still r directions in the 



equation (2.11). Therefore it is difficult to obtain the corresponding Hofstadter's butterfly 



by diagonalizing the Harper's equation (2.11): its wings are almost disappearing. 



Let me comment on the relationship between the non-commutative space and (2.11) 



Indeed the higher dimensional tight-binding model, discussed in this section, can be rep- 
resented in terms of the non-commutative torus. The coordinates of the two-dimensional 
non-commutative torus Tq are given by 

2-nid 



uv 



(2.12) 
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with e 2mS being the q-th. root of unity, 2tt9 = uj = 2-K/q. They can be written in q x q matrix 
forms, 

/o 1 \ (l \ 



u 



1 



1 



1 




V 



J2iu> 



\ 



,(q-l)iu) 



(2.13) 



Based on this description, we can consider d = 2r dimensional non-commutative torus: 

U s V t = e 2ni6st V s U t for s,t = l,---,r. (2.14) 



In general, 9 s t is a real anti-symmetric matrix. In the case of (2.11), it is simply given by a 



diagonal matrix, 2n6 s t = uj s 5 s t- This means the non-commutativity on the 2r-dimensional 
torus is introduced to each two-dimensional subspace, T| r — > x ••• x Tg Then the 
associated Hamiltonian can be written as 



T 



(2.15) 



Although, precisely speaking, we have to include a factor corresponding to the plane wave, we 
now omit these factors for simplicity. This representation means that the non-commutative 
torus operator plays a role of the translation operator with the external magnetic field even 
in the higher dimensional case. 



2.1.1 7r-flux state 

Let us comment on the 7r-flux state in higher dimensions. For the 7r-flux state [111 [T2] all 
the plaquettes have the same value, P^ u = — 1, independent of its position x, and directions 
fi, v as 

iV(z) = V^x)U v (x + fi)Ul(x + v)Ul{x) = -1. (2.16) 
Such a configuration is realized by the following Z2 C U(l) link variables, 



U^x) = rfc = (-1) 



xi H hx M _i 



(2.17) 



Then the tight-binding Hamiltonian in the second quantized form with this gauge configu- 
ration yields 

d 

X [1=1 

We now consider the free field theory for simplicity. This is almost the same as the staggered 
Dirac operator [T71[T8l[T9] 

d 1 

^staggered = ^ ^ l£)yXx {.Xx+p, ~ Xx-fi) ■ (2.19) 
x /^=1 
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Actually, by applying the transformation c x — > i xi ~^ i ~ Xd c x , the Hamiltonian (2.18) is rewrit- 
ten as 

d 

■Htight = J2J2 i ^4 {cx+fL - c x -p,) . (2.20) 

X f_l=l 



This is just the staggered fermion action (2.19), up to a constant factor. Note that the 



staggered Dirac operator is anti-Hermitian while the tight-binding Hamiltonian is Hermitian. 

There is no spinor structure in this formulation, but this staggered fermion is directly 
obtained from the naive Dirac fermion through the spin-diagonalization. See Appendix 
[A] for details. The sign factor r/ M = (— l) Xl_l ^m-i is a remnant of the gamma matrix. 
The staggered fermion enjoys an exact U(l) chiral symmetry, which is generated by e x = 
(— l) Xl ^ l_Xd . This symmetry ensures its gapless excitation. Therefore the 7r-flux state is 
also generically gapless. 

2.2 Dirac fermion models 

We then attempt to extend the Hofsdater's problem to the relativistic system. Similar 
approaches have been seen in the context of graphene [201 IZH E21 [23] . 
The generic form of the Dirac fermion action is written as 

5 D iraC = ^ \i) x+jj U^{x)P^ x ~ ^ X Ul(x)P~1p x+ p] = ^^D^yl/jy. (2.21) 

x,fj, x,y 

Here P^ stands for the spinor matrix, corresponding to a type of the lattice fermionQ The 
simplest one is given by 

^ = (2-22) 

This is just the naive lattice discretization of the relativistic fermion (25| We have other 
choices for the spinor matrix: 



• Wilson fermion [25 



1 



P± = -( 7M ±rl) (2.23) 



2 



Karsten-Wilczek fermion [271 EE] 



Pi = \lv (M = <0. P i = \{ln±ir ld ) (otherwise) (2.24) 



Borici-Creutz fermion |29|. [30] 

\{{l±ir) ltl TirT), r = i^)^ (2.25) 



For example, see [24] 



2 As well known this simple lattice discretization scheme has a problem: there exist extra massless modes 
at low energy, which are called the species doublers. To obtain a single chiral fermion we have to implement 
a much complicated scheme, e.g. domain-wall fermion, overlap fermion and so on (see, for example, a recent 
review [26]). However this problem does not concern our study because it does not affect on the spectrum of 
the Dirac fermion. 
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All these fermions include a free parameter r. We often consider the case of r = 1 for 
simplicity. Note that another lattice fermion, which is called the staggered fermion |17l [TBI 
I19j . is essentially equivalent to the naive fermion, and also the 7r-flux state as discussed in 



section 2.1 The naive and staggered fermions are transformed to each other by the spin- 



diagonalization. Thus in this paper we do not deal with the staggered fermion explicitly. 



The Dirac equation, which corresponds to the Schrodinger equation (2.8), is given by 

d 



[^(x)P+Vx +A 



Utix - fi)P~^ x - 



(2.26) 



When we write this equation (2.26) as Dip = Xip, although this Dirac operator D becomes 



non-Hermitian in general, one can define an alternative Hermitian operator, H = In 
the following we consider the spectrum of this Hermitian operator as Hip = Eip, instead of 



the original Dirac equation (2.26). 



We can also obtain the corresponding Harper's equation by taking Fourier transformation. 



In this case, the Harper's equation is slightly modified from (2.11) as 
+ ( zsin (k 2s -i - oj s Zs) (P^-i + P 2s-i) + cos (fes-l - u s z s ) (P, 



d/2 

8=1 



2s-l 



(2.27) 



Here we have 2N q energy bands due to the 7-matrix structure: they appear in a pair for the 
relativistic theory. In this case, it is again difficult to obtain a fully gapped spectrum since 
this reduced equation is not one-dimensional, but (i/2-dimensional. 



3 Non-Abelian gauge field models 

We consider another generalization of the Hofstadter problem in higher dimensions by apply- 
ing a non-Abelian gauge field as a background configuration. In this paper we concentrate 
on the four dimensional model with SU(2) gauge field. A further generalization to higher 
dimensional, and higher rank systems seems to be straightforward. 

3.1 Background configuration 

We now introduce the following SU(2) background configuration as a simple generalization 



of the U(l) background (2.2) 



Ac 



0, 



Aj = —ujjXoaj for j = 1, 2, 3, 



with aj being the Pauli matrix. Here we define filling fractions 

■El 



2ir pj 



(3.1) 



(3.2) 
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From this configuration we compute a field strength, 

F j = —UjCTj, Fjk = —LUjUkxleijk<J l ■ (3.3) 
The total background flux is given by 

J (fx TV F^F pa = J d 4 x x\ = 2ttN. (3.4) 



16vr 2 

Here this integral is taken over the four dimensional hypercubic lattice of the size L 4 , and 
we define the flux number N = mn-znz. 



Thus the link variables associated with (3.1) are given as follows, 



U (x) = 1, Uj(x) = e i0J i x ° a i = 1 coswj-xo - i&j smu)jX . (3.5) 

In this case, the translation symmetry of the lattice system yields 

x ~ x + q L CM, Xj ~ xj + 1, (3.6) 

where (/lcm is the least common multiple of (71,(72 and §3. This means that the unit cell of 
this system is also extended in only one dimension as well as the two dimensional case with 
U(l) magnetic field. 

3.2 Lattice fermion models 

Let us first study the tight-binding model with the non-Abelian gauge field configuration. 



With a Fourier basis for xq = (/lcm^Q) + the Hamiltonian (2.7) gives rise to the corre- 



sponding Harper's equation (2.11) 



3 

4>x +i + i>x -i + ( cos % cos u}jx + Oj sin/c,- sin cjjX J i/> So = E^^. (3.7) 
3=1 

This is just an effective one-dimensional model as the standard two dimensional Hofstadter 
problem. The original four dimensional model is reduced to this due to SU(2) = S 3 gauge 
symmetry. On the other hand, in this case it is also difficult to show its spectrum is totally 
gapped as the case of the U(l) background in higher dimensions. Therefore we now study 
its density of states instead of the original energy spectrum. 

Fig. [T] shows the density of states for the Harper's equation (3.7) against coo, with homo- 
geneous configuration ojj = ujq for j = 1,2,3. Although its spectrum is not totally gapped, 
one can observe a hierarchical structure of the spectrum. At the even fraction flux, e.g. 
ujo/(2tt) = 1/2, 1/4, • • •, we find a dip in the spectrum. Similar structure is discussed in the 
ordinary two dimensional Hofstadter problem at the gapless point of the filling fraction. 

In Fig. [2] we show the spectrum in the inhomogeneous flux with fixing ratios of 0/7=1,2,3- 
When flux in some directions is turned off, e.g. cJ = u;q(1,0, 0) and ujq(1, 1,0), the charac- 
teristic structure of the Hofstadter's butterfly cannot be observed in the spectrum: gap is 
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0.045 




Figure 1: Density of states for the Harper's equation (3.7) of the flux 0^=1,2,3 = ojq = 2np/q 
with g = 100, p = 0,1,2, •••99. 



completely closed. In the cases of Co = u)q(2, 1, 1) and ujq(2, 2, 1), there are some gaps in the 
spectrum and the hierarchical structure. We also find a dip at the gapless point, ujq = tt. 

Let us now comment on the relation to the non-commutative space. We introduce the 
following non-commutative torus, 



UVj = e w ^VjU for j = 1, 2, 3. 



Here U is defined as well as (2.13), while Vj=\ t 2,3 are slightly modified as 



Vj = diag ( 1, e 1 ^, • • • , e {qLCM - 1)iu1 ^ 



(3.8) 



(3.9) 



In this case the non-commutative parameter is matrix valued. Thus, the commutation 
relations between V^=i ; 2,3 cannot be written in a simple way. 

As discussed in Sec.[2j the tight-binding Hamiltonian with the external field, correspond- 



ing to the Harper's equation (3.7), can be associated with the non-commutative torus. The 



naive form of the Hamiltonian without momentum dependence is given by 

3 



U + 



v j + v j 



(3.10) 



3=1 



However, since the second part is represented as 



Vj + Vj = diag (2, 2 cos uij, • • • , 2cos(gLCM - l)^j) x 1, 



(3.11) 



we cannot involve a matrix structure in this way. Actually the SU(2) nature of the external 



field is coming through the momentum dependence in (3.7). 
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Figure 2: Density of states for the Harper's equation (3.7) with the inhomogeneous flux for 



u = 0J n: (a) n = (1, 0, 0), (b) (1, 1, 0), (c) (2, 1, 1), and (d) (2, 2, 1). 



We then investigate the Dirac fermion models as discussed in section 2.2 Applying the 



background configuration (3.5), the corresponding Dirac equation is given by 



75 



3 + 1 



+ ^2 ( ( cos fej cos UjZo + f j sin fcj sin Wj-sco) ( Pj~ — Pj 

3=1 



+ (sin fcj cos xq — (Jj cos fcj sin UjXo) I PJ + P- j j ipx 



x ■ 



(3.12) 



Here we again show the Hermitian version of the Dirac equation by multiplying the matrix 
75. When we write this Dirac operator in a matrix form, its matrix size is 8(/lcm = 2 (color) 
x 4 (spinor) x (/lcm (flux). Here color corresponds to the rank of the gauge flux of SU(2). 
On the other hand, the number of eigenvalues is given by 4<7lcM) since each spectrum is 
doubly degenerated due to the spinor structure. 



4 Summary and discussions 

In this paper we have explored some extensions of the Hofstadter problem in higher dimen- 
sions. First example is formulated with Abelian gauge configuration in d = 2r, giving rise to 
non-zero topological number. We have shown that half of the gauge potential can be triv- 
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ial by applying the Landau gauge, thus the corresponding Harper's equation is essentially 
written as r-dimensional lattice model. 

We have also pointed out that the 7r-flux state is equivalent to the staggered formalism 
of the relativistic lattice fermion. The latter is directly related to the naive Dirac fermion 
through the spin diagonalization. This means the 7r-flux state involves the chiral symmetry, 
and thus it yields massless excitation in any dimensions. 

We have then investigated SU(2) non-Abelian gauge field configuration in four dimen- 
sions. We have considered the configuration with one specific direction in four dimensions. 
In this case hopping terms for the other three directions are reduced due to SU(2) = S* 3 sym- 
metry of the gauge field. Thus we have obtained the one dimensional Harper's equation by 
utilizing the Fourier basis. We have calculated its spectrum numerically, and its hierarchical 
structure is actually observed. 

Let us now comment on possibilities of future works along this direction. The two 
dimensional Hofstadter problem is essentially related to the quantum group [3~T| l32l [33], [20] : 
the Harper's equation is directly regarded as the Baxter's equation for the one dimensional 
U q {s\2) model. Thus it is interesting to explore the corresponding quantum group structure to 
the generalized Hofstadter problems discussed in this paper. In particular, the non-Abelian 
version of Harper's equation includes the matrix-valued coefficient. This corresponds to q- 
parameter in the two dimensional case, thus it is natural to investigate a kind of quantum 
group with matrix-valued ^-parameter. 

Next is the lattice study with various kinds of lattice fermions, i.e. Wilson [25], staggered 
dZldiKin], staggered- Wilson O E3 E3 EH] , minimal-doubling (23 |2Sl HE E3 ESI 113 E] , 
domain- wall |42^ H3] 133] and overlap fermions [15] . They were originally introduced to tackle 
the difficulty of the chiral fermion on the lattice, but they themselves are interesting as 
statistical lattice models: some of them are actually investigated in the context of condensed- 
matter physics, for example, graphene, 7r-flux state, topological insulator /superconductor 
and so on. Thus we hope the Hofstadter problem formulated with these lattice fermions are 
relevant to realistic condensed-matter physics. 
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A Spin diagonalization 

We now show that there is an alternative expression of the d-dimensional naive Dirac fermion 
without spinor matrix structure. It is given by diagonalizing the corresponding 7 matrices. 
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Let us start with the naive Dirac fermion on the lattice, 

1 



5 = E 



X,fl 



(A.l) 



Then, introducing the field Xx defined as 

/ Xj Xo Xl 7 - X\ Xo x, 

Vx = l d d • • -72 7l Xx, Vx = Xxlx ll ""Id 



(A.2) 



we can represent the naive lattice fermion (|A.1|) in the following form, 
1 



5 = E 



X,/l 



:V»( x )Xx (Xx+n ~ Xx-p,) + mxxXx 



with nn(x) = (-1) E "<"^. (A.3) 



Remark there is no spinor structure in this expression. In other words, the spinor matrix 
is diagonalized in this basis. This means the naive Dirac fermion can be rewritten in terms 
of one-component fermionic field. This is just the staggered formalism of the relativistic 
fermion [THIIIS]- 
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